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Abstract 

ON 

We update a previously-proposed set of supersymmetric benchmark scenarios, taking 
into account the precise constraints on the cold dark matter density obtained by combining 
WMAP and other cosmological data, as well as the LEP and b — > s'j constraints. We assume 
that R parity is conserved and work within the constrained MSSM (CMSSM) with universal 
soft supersymmetry-breaking scalar and gaugino masses m and m^. In most cases, the relic 
density calculated for the previous benchmarks may be brought within the WMAP range by 
reducing slightly mo, but in two cases more substantial changes in mo and mi/2 are made. 
Since the WMAP constraint reduces the effective dimensionality of the CMSSM parameter 
space, one may study phenomenology along 'WMAP lines' in the (mi/2,mo) plane that have 
acceptable amounts of dark matter. We discuss the production, decays and detectability of 
sparticles along these lines, at the LHC and at linear e + e~ colliders in the sub- and multi-TeV 
ranges, stressing the complementarity of hadron and lepton colliders, and with particular 
emphasis on the neutralino sector. Finally, we preview the accuracy with which one might be 
able to predict the density of supersymmetric cold dark matter using collider measurements. 
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1 Introduction 



One of the crucial topics in the planning of analyses of data from experiments at present 
and future colliders is the search for supersymmetry [1]. Even the minimal supersymmetric 
extension of the Standard Model (MSSM), which conserves R parity, has over 100 free 
parameters once arbitrary soft supersymmetry-breaking parameters are allowed. For this 
reason, much attention is focussed on the constrained MSSM (CMSSM), in which the soft 
supersymmetry-breaking scalar masses m , gaugino masses mi/ 2 and trilinear parameters A 
are each assumed to be universal at some high input scale, as in minimal supergravity and 
other models. These CMSSM parameters are constrained principally by the absence at LEP 
of new particles with masses ^ 100 GeV, by the agreement of b — > 57 decay with Standard 
Model predictions, by the measurement of the anomalous magnetic moment of the muon, 
— 2, and by the range allowed for relic cold dark matter, Q C dm- 

Benchmark supersymmetric scenarios have a venerable history [2]. A couple of years 
ago, a new set of benchmark supersymmetric models was proposed in [3], consistent with all 
the above experimental constraints, as well as the cosmological constraint on VLcdm- Sub- 
sequently, there has not been any significant change in the LEP limits on supersymmetric 
particles [4, 5] , we are still looking forward to more sensitive searches at the Fermilab Teva- 
tron collider, the situation of b 57 decay has changed little [6], and the interpretation of 
the 5^ — 2 measurements remains unclear [7] . Various other benchmark scenarios have been 
proposed, in particular some supplementary points and lines in the CMSSM and in models 
with different mechanisms for supersymmetry breaking [8] , benchmarks for supersymmetric 
Higgs physics [9], and scenarios in which the prospects for the Fermilab Tevatron collider 
are more favourable [10, 11]. 

The road(w)map for supersymmetric phenomenology has, however, been altered signifi- 
cantly by the recent improved determination of the allowable range of the cold dark matter 
density obtained by combining WMAP and other cosmological data: 0.094 < Qqdm < 0.129 
at the 2-0" level [12]. This range is consistent with earlier indications, but more precise. 
Within the MSSM with conserved R parity, if one assumes that most of the cold dark mat- 
ter consists of the lightest supersymmetric particle (LSP) and identifies this with the lightest 
neutralino \ [13], this WMAP constraint reduces the dimensionality of the parameter space. 
In the CMSSM, in particular, whereas generic regions of the (m^jfno) planes would pre- 
viously have been allowed for fixed values of tan/9 and Aq [14], now only thin strips are 
permitted [15, 16] 1 . 

1 Some authors had also considered an analogous narrow range for Q, x h 2 before this was mandated by the 
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It is of course possible to relax the restrictions to these narrow 'WMAP lines' in various 
ways. One could consider models that violate R parity, though these would also be subject 
to (different) astrophysical and cosmological constraints. Or one could consider alternatives 
to gravity-mediated supersymmetry breaking, such as anomaly- or gauge- mediated super- 
symmetry breaking, in which the sparticle spectra could be radically different, and/or the 
LSP might be a gravitino [8] . However, we do not discuss such models in this paper. 

Here, we consider first the minor modifications of the previous CMSSM benchmark sce- 
narios that bring them on to the WMAP lines. In most cases, this requires only a small 
change in m , keeping the same as before, as seen in Table 1. An exception is the 
previous benchmark H, which had been chosen at the end of one of the coannihilation 'tails', 
at the largest possible value of m^. The WMAP constraint not only reduces the possible 
spread in m for fixed my 21 but also reduces the allowable range of mi/ 2 . Therefore, our 
modification of benchmark H has reduced values of both m 1 / 2 and m . Another example is 
benchmark M, which has also been shifted significantly because the rapid-annihilation 'fun- 
nel' for tan/5 = 50 and /i > is affected substantially by the WMAP restriction on Vt x h 2 . 
More detailed discussions of these updated benchmark points are given in Section 2, where 
we also discuss the (generally small) extent to which the decay branching ratios differ from 
the previous versions of the benchmark points 2 . 

Then, in Section 3, we discuss systematically how the phenomenology of supersymmetric 
models varies with along the WMAP lines. We first present parametrizations of these 
lines appropriate for the SSARD [18] and ISASUGRA 7.67 [20] codes. We then discuss the 
patterns of decays of various selected sparticles. We show in particular that X2 t\t decays 
are generally important, though other % 2 — > ££ decays (where £ = e, /i) are also important 
for low mi/2 at low tan (3. The decays X2 — > X^, X% are never very large on the WMAP lines, 
despite being kinematically allowed for most values of my 2 - We then discuss the average 
numbers of Z, h and r particles produced per sparticle production event at the LHC. This 
study confirms the importance of the r signature, with trilepton signatures also looking 
promising. We discuss the extent to which our various benchmark scenarios sample generic 
features along the WMAP lines. 

Section 4 contains a discussion of sparticle detectability at various accelerators, starting 
with the LHC 3 . We find that the lightest MSSM Higgs boson and all the squarks are in prin- 
WMAP data [17]. 

2 Where it seems necessary to avoid confusion with the previous versions of these benchmark scenarios [3], 
we denote the updated versions with primes: A', B', etc.: otherwise we retain the same notation. 

3 Wc have nothing to add here to the discussion in [3] of the prospects for detecting supersymmetry at 
the Fermilab Tevatron collider, referring the reader to [10, 11] for alternative discussions. 
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Supersymmetric spectra in post-WMAP benchmark scenarios 
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Table 1: 

Proposed post-WMAP CMSSM benchmark points and mass spectra (in GeV), as calculated 
using SSARD [18] and FeynHiggs [19], using the one-loop corrected effective potential com- 
puted at the electroweak scale and one-loop corrections to the chargino and neutralino masses. 
We recall (in parentheses) the values of mi/ 2 and m used in [3], in cases where they differ. 
As in [3], exact gauge coupling unification is enforced, resulting in the predictions for a s {mz) 
that are shown in units of 0.001. We use A = 0, m^m;,)^ = 4.25 GeV and mt = 175 GeV 
for most of the points, with the exceptions of points E' and F', where m t = 171 GeV is used. 
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ciple observable along the complete WMAP lines, except along parts of the rapid-annihilation 
'funnels', whereas sleptons and other neutralinos and charginos are detectable only in the 
lower parts of the WMAP ranges of mi/ 2 . This discussion is extended in Subsection 4.2 to 
linear e + e~ colliders with E C m = 0.5, 1.0 [21] and 3.0 or 5.0 TeV [22]. Confirming previous 
studies [21, 3, 8], we see that a machine with E C m = 0.5 TeV would be able to explore 
supersymmetry in the lower parts of the WMAP ranges of mi/ 2 , whereas a machine with 
Ecm = 1-0 TeV would be able to explore supersymmetry along (almost) the entire WMAP 
lines. CLIC with E C m = 3.0 TeV would be able to complete the spectrum of electroweakly- 
interacting sparticles along the entire WMAP lines, as well as furnish detailed measurements 
of squarks for a large fraction of the allowed range of mi/ 2 [23]. We also discuss in more detail 
the observability of neutralinos at the different colliders surveyed, with particular emphasis 
on CLIC studies. 

Finally, in Section 5 we review some conclusions from our results, addressing in particular 
the question whether collider measurements will have the potential to determine if the LSP 
constitutes most of the cold dark matter in the Universe [13, 24]. 

2 Updated CMSSM Benchmark Points 

2.1 Improved Choices of Supersymmetry-Breaking Parameters 

Since the laboratory constraints on the CMSSM have not yet changed substantially since 
the termination of LEP [4] , the only reason for updating the benchmark points we proposed 
previously [3] is the refined estimate of VlcDuh 2 provided by combining the new WMAP 
data with those previously available [12]. Assuming that most of the cold dark matter 
is composed of LSPs \i previously we allowed 0.1 < Q x h 2 < 0.3 (rather conservatively), 
whereas the WMAP analysis now allows only the range 

0.094 < tt x h 2 < 0.129 (1) 

at the 2-a level [12]. We see in Table 2 of [3] that all the previous benchmark points 
yielded relic densities Q x h 2 above the range (1). Since the relic density calculations have not 
changed significantly since [3] in the regions of interest, the previous benchmark points must 
be moved. 

Of the 13 benchmark points proposed previously, 5 (B, C, G, I, L) were in the 'bulk' 
regions at low mi/ 2 and m , 4 (A, D, H, J) were along the coannihilation 'tails' extending 
to larger mi/ 2 [25], 2 (K, M) were along rapid-annihilation 'funnels' where both mi/2 and 
m may grow large [26], and 2 (E, F) were in the 'focus-point' region at very large m [27]. 
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The WMAP constraint (1) has slimmed the 'bulk' region down considerably, but only minor 
reductions in m are needed to reduce the relic density of points (B, C, G, I, L) into the range 
(1), as shown in Table 1. The coannihilation 'tails' are also much thinner than they were 
before, and the required reductions in m for points (A, D, J) are also small. The previous 
point H is an exception to this rule, since it was chosen at the extreme tip of a coannihilation 
'tail'. In this case, since the upper limit in (1) is considerably reduced compared with that 
assumed in [3], the allowed upper limit on mi/ 2 has been reduced substantially [15]. Hence 
we must make reductions in both mi/ 2 and m for point H, as also shown in Table 1. In the 
case of point K, the rapid-annihilation 'funnel' has become thinner following WMAP, and a 
minor adjustment in mi/ 2 is sufficient, but in the case of point M, the 'funnel' has changed 
substantially, and both m^j and m had to be changed. 

The two focus-point benchmarks (E, F) may also be adapted to the WMAP constraint (1) 
with small changes, as also shown in Table 1, which is based on the SSARD code. However, 
we take this opportunity to underline the extreme delicacy of model calculations in this 
region, as may be seen by comparing the results of different codes [18, 20, 28, 29] in their 
successive releases 4 . The various codes do agree that there is a strip in the 'focus-point' 
region where tt x h 2 is in the WMAP range (1), but the error in its location due to the 
experimental uncertainty in the mass of the top quark m t (for example) is much larger 
than the intrinsic width of the WMAP strip in this region. We have chosen a pole mass 
m t = 171 GeV for benchmarks E and F, as opposed to the choice m t = 175 GeV made for the 
other benchmarks, which would have required much larger values of m in this focus-point 
region, namely m = 2550, 5030 GeV for the SSARD code, respectively. The new DO value 
of m t ~ 179 GeV [30], if confirmed, would push the 'focus-point' region up to still larger 
m , whatever code is used: specifically to m = 5800, 9070 GeV for benchmarks E and F, 
respectively, if the SSARD code is used 5 

There are also significant theoretical uncertainties associated with higher-order effects [31], 
that are reflected in differences between different codes in their various versions. For exam- 
ple, for m t = 175 GeV at benchmarks E and F (with mi/ 2 = 300, 1000 GeV, respectively), 
the ISASUGRA 7.51, [ISASUGRA 7.67], [SUSPECT 2.10], [SUSPECT 2.11] codes require m = 
1530,3450 GeV, [m = 3590,6260 GeV], [m = 2350,4110 GeV], [m = 2590,3850 GeV], 
respectively. We see that the predictions of successive versions of the SUSPECT code [28] 
do not vary much, and agree quite well with SSARD [18], whereas there has been significant 
evolution in the predictions of ISASUGRA [20]. 

4 See also the discussion in [3]. 

5 This would put all sfermions beyond the reach of all the colliders discussed here. 
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In view of these uncertainties in the focus-point region, we concentrate in the following 
on the remaining updated benchmark points proposed in Table 1. These are located on the 
'WMAP lines' shown in Fig. 1(a) for \l > and tan (5 = 5, 10, 20, 35 and 50, and (b) for \l < 
and tan/? = 10 and 35 6 . We recall that, for given values of tan (3, mi/ 2 and the sign of /x, 
lower values of m generally have values of Q x h 2 below the WMAP range 7 , whereas higher 
values of m generally have values of fi x h 2 that are too high, and are therefore unacceptable 
unless one goes beyond the CMSSM framework used here. At the ends of the WMAP lines', 
smaller values of mi/ 2 are excluded by either m h [5, 19] and/or b — > [6], and larger values 
of mi/2 have values of £l x h 2 above the WMAP range and/or the LSP is the charged f\. We 
note that, as was commented in [3], most of the proposed benchmark points for // > yield 
a value of — 2 that lies within 2 a of the present experimental value based on e + e~ data, 
corresponding to the lighter (pink) regions of the strips in Fig. 1(a). However, we do not 
impose this as a requirement on the benchmark points, as exemplified in Fig. 1(b) for \i < 0. 

We see also in Fig. 1 that the majority of the benchmark points lie on the portions of 
the WMAP lines with lower values of m 1 /2, making them more interesting for the early 
stages of LHC operation, and for any sub-TeV linear e + e~ collider. Generally speaking, less 
fine-tuning of parameters is required in this region than in the focus-point, funnel and tail 
regions [32]: see also the discussion in [3]. 

2.2 Discussion of Spectra 

In order to study the detectability of MSSM particles and the possible accuracies in measure- 
ments for the proposed benchmarks, it is essential to relate the input CMSSM parameters, 
defined in Table 1 for the SSARD program, to those necessary to obtain an equivalent MSSM 
spectrum in Monte Carlo generators 8 . Accordingly, we have tuned the inputs for ISASUGRA 
7.67 to reproduce as accurately as possible the spectra given in Table 1, by scanning m 
and tan/5 to identify the parameter pairs minimizing the sum of the relative differences of 
particle masses |m SSARD — m ISASUSY |/M SSARD . In view of their importances for the calculation of 
fl x h 2 , we have treated specially the mass splitting jm^ — m x |, which is important along the 
coannihilation 'tail', and also jm^ — 2m x |, where it becomes important in the 'funnel' regions. 
We have given each of them 30 % of the total weighting when optimizing the Isasugra 7 . 67 
parameters in the relevant (mi/ 2 ,mo) regions. The results of the optimization are given in 

6 There would be corresponding WMAP lines in the focus-point regions, whose location is uncertain, but 
which would generally lie at larger values of too. 

7 Therefore, in this region the LSP might constitute only a part of the cold dark matter. 
8 See also the discussion in [3]. 
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Figure 1: The shaded strips display the regions of the (mi/2, mo) plane that are compatible [15] 
with 0.094 < VL x h 2 < 0.129 in the 'bulk', coannihilation 'tail', and rapid- annihilation 'funnel' 
regions, as well as the laboratory constraints, for (a) \i > and tan j3 = 5,10,20,35 and 
50, and (b) for fi < and tan/3 = 10 and 35. The parts of these 'WMAP lines' for \x > 
compatible with — 2 at the 2-a level have lighter (pink) shading [7]. The updated post- 
WMAP benchmark scenarios are marked in red. Points (E',F') in the focus-point region 
have larger values of mo. 
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Table 2. 

In general, good matching between the predictions of the two codes can be found with 
only moderate shifts of the input parameters, and we have chosen only to adjust m in most 
cases. The typical average relative difference in the benchmark point masses is of the order 
of a few percent. However, occasionally individual masses may exhibit discrepancies of up to 
15%, and ensuring compatible mass splittings AM requires a systematic shift in m to higher 
values, particularly at large tan (3, as seen in Fig. 2, where the cases tan f3 — 5, 10, 20, 35 and 
50 are exhibited for /i > 0, and tan/5 = 10 for /i > 0. In addition to increasing m , here 
we find a better correspondence between SSARD and ISASUGRA 7.67 if tan (3 is allowed to be 
reduced, as done for benchmarks (K, L, M) in Table 2. 

We have further studied whether the resulting ISASUGRA 7 . 67 versions of the benchmark 
points yield values of VL x h 2 that are indeed compatible with the WMAP data. The first line 
of Table 3 shows the values of VL x h 2 calculated using the SSARD [18] code at the updated 
benchmark points shown in Table 1, which are well within the WMAP range. The second 
line shows the corresponding values of fi x h 2 calculated using the Micromegas code [33] 
interfaced with the ISASUGRA 7.67 code that is used to calculate the spectra in Table 2. 
The agreement between SSARD and Micromegas/ISASUGRA 7.67 is generally good in the 
bulk and coannihilation tail regions, with the latter falling within the WMAP range, except 
for point H (which is right at the tip of the coannihilation tail) and point J (where the 
discrepancy is only marginal). The Micromegas/ISASUGRA 7.67 code gives results that are 
outside the WMAP range in the focus-point and at the upper edge of the funnel region 
(points E, F and M), reflecting the fact that the calculation of Vl x h 2 is notoriously delicate 
in these regions [32] and that the tuning of the mass spectra between the two different codes 
becomes difficult. 

For completeness, we also present in Table 2 the values of the b — > 57 branching ratio and 
the supersymmetric contribution to — 2 calculated using SSARD. These may be compared 
with any evolution in the future experimental values of these quantities. As already men- 
tioned, the situation with regard to any possible non-Standard Model contribution to the 
experimental value of — 2 is sufficiently volatile that we do not use this as a criterion for 
selecting benchmarks, though it is used in Fig. 3 to order the benchmark points, as described 
below. 
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Supersymmetric spectra in post-WMAP benchmarks 
calculated with ISASUGRA 7.67 
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T2 


425 


191 


291 


376 


1534 


3485 


290 


674 


312 


579 


1251 


420 


1504 


V T 


411 


167 


273 


360 


1532 


3478 


266 


657 


278 


558 


1239 


387 


1492 


u L , c l 


1248 


558 


859 


1103 


1639 


3923 


814 


1885 


778 


1554 


2722 


1001 


3670 


Ur, Cr 


1202 


542 


830 


1064 


1637 


3897 


787 


1812 


754 


1497 


2627 


969 


3528 


dL, s L 


1250 


564 


863 


1107 


1641 


3924 


818 


1887 


783 


1556 


2723 


1004 


3671 


dR, Sr 


1197 


541 


828 


1059 


1638 


3894 


786 


1804 


752 


1491 


2615 


967 


3509 


tl 


958 


411 


653 


860 


1052 


2647 


617 


1477 


584 


1207 


2095 


753 


2857 


t 2 


1184 


576 


837 


1048 


1387 


3373 


792 


1753 


748 


1428 


2366 


920 


3231 




1147 


514 


789 


1015 


1375 


3356 


737 


1719 


677 


1377 


2297 


844 


3149 


b 2 


1181 


535 


816 


1043 


1602 


3816 


770 


1761 


725 


1423 


2349 


904 


3217 



Table 2: Proposed post-WMAP CMSSM benchmark points and mass spectra (in GeV), as 
calculated using ISASUGRA 7.67 [20] and adapting the values of mo and tan/3 (when it is 
large) to give the best fit to the SSARD spectra shown in Table 1, as described in the text. 
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Figure 2: Comparison of the mass differences AM = |m fl — m x | a/ong the WMAP lines, as 
calculated in the SSARD code (dots) and ISASUGRA 7.67 using the same values ofm Q (dashed 
lines), as functions of for /x > and tan/5 = 5,10,20,35,50 and /x < 0, tan/3 = 10. 
W'e see iTte need to s/ii/t m systematically, particularly at larger m\/2 and tan/3. This has 
been done for the solid lines, which are improved ISASUGRA 7.67 fits to the WMAP lines 
described in the text. 
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Relic density Vt x h 2 ,b — > S7 and — 2 in post-WMAP benchmark scenarios 
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4.5 
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Table 3: Comparison of Q x h 2 for the benchmark points in Table 1 computed with the SSARD 
code [18] (top line), and the Micromegas code [33] interfaced with ISASUGRA 7.67 [20] (sec- 
ond line) using the fitted parameters shown in Table 2. The third and fourth lines show the 
values of b — > S7 and g^ — 2 calculated using SSARD. 



2.3 Observability at Different Accelerators 

As an important ingredient in the subsequent analyses of sparticle observability at different 
colliders, we first report a comparison of the principal sparticle branching ratios at the up- 
dated benchmark points with those found previously [3] at the original benchmark points. 
The new branching ratios have been calculated using ISASUGRA 7 . 67, whereas previously 
they were calculated using ISASUGRA 7.51 [20], and we comment on the differences asso- 
ciated with the improvement in ISASUGRA, as opposed to the changes in the benchmark 
parameter choices 9 . Since the differences in branching ratios are not large, in general, we 
limit ourselves to commenting on instances where the more significant differences occur. 

• In the case of point A, the most notable difference is an increase in B(b 2 — > b\) from 39 
to 62 %, with corresponding decreases in B(b 2 — > t + x~ , ti + W~). There is also an increase 
of Xa ~ * from 23 to 29 % at the expense of x* ~ *■ Wxf- These changes are not related 
to the new parameter choice, but to the differences in the ISASUGRA versions. 

• In the case of point B, there is again an increase in B(b 2 — > 6%), from 15 to 25 %, and 

a decrease in B(b 2 — > t\ + W~). There is also an increase in B(e~L —> &X) from 47 to 87 %, 

accompanied by decreases in B{e^ — > ex2, v + xf )■ We also note that the B{pi — > l~x + ) now 

vanish, so that the i>i are invisible. There is also an increase in B(f 2 — > r%) from 51 to 84 %, 

which decreases B{f 2 — > v T + x~, r + X2)- Finally, there are decreases in B(x 2 — > fi + r) from 

84 to 42 %, which increases the branching ratio for the invisible channel B{x2 v + and 

in Bi^x^ 1 t~\ + v t) from 96 to 36 %, which increases B{x ± + The changes for the 

sleptons are mostly because they are lighter at the updated benchmark point with smaller 

mo, and those for b 2 are related to the differences in the ISASUGRA versions. Both effects 

9 For the purposes of our discussion, one of the most important upgrades to ISASUGRA has been the 
inclusion of one-loop radiative corrections to sparticle masses. 
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influence the gaugino branching ratios. 

• At benchmarks C and D, B(b 2 — > bx) increases from 13 to 24 % and from 20 to 44 %, 
respectively, with corresponding decreases in B(b 2 — > t\W~) and B(b 2 — > t + x 1 * 1 )- At point 
C, decreases are observed of %2 — ► T\T from 23 to 13 % and of — > ^z/ from 21 to 11 %, 
compensated by increases of X2 — > vv and xt ~^ vl, respectively. All changes are mostly 
related to differences in the ISASUGRA versions. 

• At the updated points E and F, the most notable effect is an increase of the branching 
ratios of squarks to gluinos, due to the increase in m and hence of the squark masses at the 
updated benchmark points. For instance, the B{b\ — > bg) changes from 45 to 52 % at E and 
from 22 to 38 % at F. They are compensated by a reduction of the decays to charginos and 
neutralinos. Moreover, at point E the decays of sleptons into xt tend to decrease and those 
into xt to increase correspondingly. 

• At the updated benchmark G, there are increases in B(h -> txt), from 13 to 22 %, 
with accompanying decreases in other modes, mainly in B(b\ — > t + xti b + X2>- A decrease 
is observed of X2 — > T\T from 82 to 62 % and of xt ~* from 81 to 57 %, compensated 
by an increase of X2 vv and xt respectively. The changes are mostly due to the 
differences in the ISASUGRA versions. 

• Although the updated point H has a significantly lighter spectrum than previously, no 
major changes (> 5-6 %) are observed in the decay branching ratios. 

• For the old point I the decay bi — > txt was kinematically forbidden. For the updated 
point it is allowed and has a branching ratio of 11 %. Also, B{u e — > vx) increases from 
61 to 73 % and B{v e — > ext) decreases correspondingly, mostly due to the change in the 
ISASUGRA version used. 

• At point J the main changes are a decrease of B(x2 — > Tit) from 82 to 66 % and of 
B(xt T\v) from 82 to 64 %, mostly compensated by an increase of B(x2 — > vv) and 
B(xt — >■ vl), respectively. The changes are mainly due to the differences in the ISASUGRA 
versions. 

• At the updated benchmark K, there is an increased branching ratio for the direct 
decay of cl and v e to x from 26 to 35 % and 27 to 37 % respectively, at the expense of 
the decays to xt- The increase of B{v T — > Wt\) from 28 to 40 % is accompanied by a 
decrease of B(y T — > rxt)- Also, the decays of X2 and xt to ?i, which were previously 
forbidden kinematically, are now allowed by the new parameter values and reach 43 and 41 
%, respectively. A reduction of B(xi Wxt) from 56 to 47 % is also observed, associated 
with the differences in the ISASUGRA versions. 

• The updated point L is characterized by a decrease of the branching ratios of u e and v T 
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to the Xi, compensated by an increase of B(v e — > vy) from 35 to 43 % and of B{v T — > W^Ti) 
from 59 to 77 %. The decay branching ratios of X3 and %4 to t\t are increased from 19 to 27 
% and 13 to 20 %, respectively, with a corresponding reduction of the decays to Xi ■ These 
changes are associated with the updated CMSSM parameter values. 

• At the updated point M, the main change is a smaller value of m , which affects many 
branching fractions. We only mention some of the major changes. The direct decay to x is 
increased for e L and V e from 34 to 74 % and from 34 to 77 %, repectively, at the expense of 
a reduced decay to the chargino and the heavier neutralino. The B(x2 — ► h + x) decreases 
from 45 to 12 %, whereas B(x2 — > T\ + r) increases, and B(xf — > W ± + x) decreases from 
48 to 14 %, whereas B(xf — > n + i/ T ) increases. 

The branching ratios at the updated benchmark points not mentioned above are not 
significantly different from those at the original versions of the points. Among all the above 
effects, perhaps the most significant is the invisibility of the i>e at the updated version of 
point B. 

Combining the above information, we now present in Fig. 3 an updated comparison of 
the numbers of different MSSM particles that should be observable at different accelerators 
in the various benchmark scenarios [3], ordered by their consistency with — 2 as calculated 
using e + e~ data for the Standard Model contribution [7] . We re-emphasize that the qualities 
of the prospective sparticle observations at hadron colliders and linear e + e~ colliders are of- 
ten very different, with the latters' clean experimental environments providing prospects for 
measurements with better precision. Nevertheless, Fig. 3 already restates the clear message 
that hadron colliders and linear e + e~ colliders are largely complementary in the classes of 
particles that they can see, with the former offering good prospects for strongly-interacting 
sparticles such as squarks and gluinos, and the latter excelling for weakly-interacting sparti- 
cles such as charginos, neutralinos and sleptons. We discuss later the detailed criteria used 
for assessing the detectabilities of different particles at different colliders. 

3 WMAP Lines 

As has already been mentioned, in view of the reduction in dimensionality of the CMSSM 
parameter space enforced by WMAP [15], one may progress beyond the previous approach 
of sampling, more or less sparsely, the CMSSM parameter space. In particular, as proposed 
in [8] but with a different attitude towards the cosmological density constraint, one may 
explore, more or less systematically, the CMSSM phenomenology along the lines in parameter 
space shown in Fig. 1. Again as proposed in [8], one may then focus special attention 



13 



gluino ^™ squarks ^™ sleptons ^™ X 

Post-WMAP Benchmarks 



H 



73 30 



& 20 r 

^ 10 

3 o 

03 

O 30 
® 20 



LHC 




I | | | 



.1 



LBG I C J HMAE FKD 







LC 1.0 TeV 



F l ■ I ■ 



LBG I C J HMAE FKD 




LBG I C JHMAE FKD 



30 
20 
10 




LHC+LC 1 TeV 




LBG I C JHMAE FKD 



30 
20 
10 





LBG I C J HMAE FKD 




LBG I C JHMAE FKD 



Figure 3: Summary of the numbers of MSSM particles that may be detectable at various 
accelerators in the updated benchmark scenarios. As in [3], we see that the capabilities of the 
LHC and of linear e + e~ colliders are largely complementary. We re- emphasize that mass 
and coupling measurements at e + e~ colliders are usually much cleaner and more precise 
than at hadron-hadron colliders such as the LHC, where, for example, it is not known how 
to distinguish the light squark flavours. 
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on particular points along these lines. The updated benchmark points discussed in the 
previous section are examples and, as we discuss below, they manifest some of the distinctive 
possibilities that may appear along the WMAP lines. However, for certain purposes it may 
be helpful to select a few additional points that exemplify other generic possibilities, as we 
discuss later. 

3.1 Parametrizations of WMAP Lines 

As seen in Fig. 1, the WMAP lines have a relatively simple form for tan /3 ^ 35. It is there- 
fore possible to provide simple quadratic parametrizations for them that are a convenient 
summaries of the correlations between m^j and mo that are enforced by WMAP within the 
CMSSM framework adopted here. Using the SSARD code, we find the following convenient 
parametrizations for the fi > lines: 

= -5.46 + 0.206m 1/2 + 4.02 x 10' 6 m 2 1/2 , (475 < m 1/2 < 960) (2) 

= 11.97 + 0.163m V2 + 4.90 x l0' 5 m 2 1/2 , (300 < m 1/2 < 940) (3) 

= 50.89 + 0.144m V2 + 6.77 x 10- 5 m 2 1/2 , (270 < m 1/2 < 940), (4) 

= 90.24 + 0.219m 1/2 + 5.70 x 10" 5 m? /2 , (315 < m 1/2 < 1070) (5) 

where the masses are expressed in units of GeV, and the numbers in parentheses specify 
the allowed ranges of mi/ 2 . When tan/3 = 50, a rapid-annihilation funnel appears, and 
parametrizing the regions allowed by WMAP becomes more complicated. In the case of 
tan (3 = 50, one may conveniently use: 

m = 140.6 + 0.33m 1/2 + 6.3 x 10~ 5 m 2 1/2 (350 < m 1/2 < 1550), 

m = -1941 + 1.65m 1/2 (1700 < m 1/2 < 1800), (6) 

where the masses are again expressed in GeV units. For jj, < 0, we propose the following 
convenient parametrizations: 

tan/3=10: m = 15.0 + 0.1Qm 1/2 + 4.2 x 10 -5 m^ 2 , (400 < m 1/2 < 940) (7) 
tan/3 = 35: m = 39.6 + 0.76m 1/2 - 1.9 x 10" 5 m? /2 , (790 < m 1/2 < 1650), (8) 

where for tan /3 = 35 we parametrize only the top branch allowed by WMAP in the (mi/ 2 , m Q ) 
plane shown in panel (b) of Fig. 1, which includes the updated point K. 

As mentioned earlier when discussing the specific benchmark points, the appropriate 
values of m as functions of mi/ 2 must be re-evaluated for the ISASUGRA 7.67 code. A 
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similar minimization procedure to that discussed earlier has been repeated, parametrizing 
the shifts w.r.t. the values obtained with the SSARD code by a linear function of m^, which 
have then been added to (5, 6, 8). The resulting parametrizations are: 

= -9.96 + 0.197m 1/2 + 4.02 x l(T 6 m? /2 , (475 < m 1/2 < 960) (9) 
= 13.57 + 0.142mi /2 + 4.90 x l^m\ /2 , (300 < m 1/2 < 940) (10) 
= 50.39 + 0.142m 1/2 + 6.77 x l(}- 5 m 2 1/2 , (270 < m 1/2 < 940), (11) 
= 89.56 + 0.239m 1/2 + 5.70 x 10~ 5 m 2 1/2 , (300 < m 1/2 < 1050). (12) 

These have been used to produce the solid lines in Fig. 2. 

We note in passing that, as the mass difference AM = \m fl — m x \ — > towards the 
high-mi/ 2 ends of the WMAP lines, these have portions at high mi/ 2 where AM < m T . In 
these portions, the fi NLSP decays via a virtual r, resulting in the dominance of four-body 
decay modes such as fi — > vxtv and v\qq. In this case, the f\ is stable on the scale of 
the size of the detector. This observation has important implications for f\ observability at 
different colliders, as we discuss below. 
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3.2 Discussion of \2 Decay Branching Ratios 

One of the key particles appearing in sparticle decay chains is the second neutralino X2, 
whose branching ratios are quite model-dependent and have significant impact on sparticle 
detectability at future colliders. Moreover, \2 decays play crucial roles in reconstructing 
sparticle masses via cascade decays. Therefore, we now use ISASUGRA 7.67 to discuss how 
the principal branching ratios of the \2 vary along the WMAP lines, noting several significant 
features that are important for phenomenology. 

As seen in Fig 4(a) for tan/3 = 5 and /i > 0, \i decays into fr and other ££ modes 
dominate among the decay modes of interest. Both tyt and r 2 r contribute, the latter in- 
creasing with mi/2 while the former decreases. The dip in the the ££ branching ratio when 
mi/2 ~ 300 GeV reflects a similar switch between the £r£ modes at low mi/ 2 and the £ L £ 
modes at large mi/ 2 10 . The most important other decay modes are the invisible X2 — > vv 
decays. The decay modes \2 — * X% an d X2 — > xh exhibit clear thresholds at mi/ 2 ~ 270 GeV 
and 310 GeV, respectively, reflecting the opening up of these two channels. We note that 
benchmark A is in the region of m^j where the branching ratio for X2 — > X^ 1 is already 
declining from its peak value ~ 6 %, while that for % 2 — > \Z has already been reduced from 
its peak value ~ 0.6 %. 

10 We note in passing that the branching ratio for i — > \(. exceeds 50 % for all the WMAP range of m^. 
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Figure 4: Dominant branching ratios of the next-to-lightest neutralino X2 as functions ofmi/2 
along the WMAP lines for (a) tan/? = 5, (b) tan (3 = 10, (c) tan/? = 20, (d) tan (3 = 35 
and (e) tan/3 = 50, all for \i > 0, and (f) tan/3 = 10 for \x < 0. The locations of updated 
benchmark points along these WMAP lines are indicated, as are the upper limits on m 1( / 2 in 
panels (a, b, c) and (f). 
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Similar features are exhibited in Fig 4(b) for tan/3 = 10 and fj, > 0, with the peak in the 
branching ratio for X2 — * xh reduced to ~ 5 % and that for \i ~^ X% increased to ~ 0.9 
%. The updated benchmark point B with mi/ 2 = 250 GeV is below both thresholds, whilst 
benchmark C has a near-maximal branching ratio for X2 - ► xh- 

Panel (c) of Fig. 4 shows the corresponding \2 decay branching ratios for tan (3 = 20 
and fx > 0. We see that the fr mode is larger than for tan (3 = 10, reflecting the larger r 
Yukawa coupling. The dip in the 11 has moved to larger toi/ 2 ~ 400 GeV. The thresholds in 
X2 X% and X2 —> xh are again prominent, with the former branching ratio now reaching 
slightly more than 1 %. Benchmark G is located near the peak in X2 — > xh and the dip in 
X2 — > it-, whilst benchmark H has almost equal branching ratios for 11 and fr, and smaller 
branching ratios for xh and particularly \Z . 

When tan (3 = 35 and // > 0, shown in panel (d) of Fig. 4, \2 — * tt becomes the dominant 
branching ratio for all values of mi/ 2 allowed by WMAP. Apart from this, the most noticeable 
feature is the relative suppression of the X2 — ► M branching ratio, whose dip has now moved 
up to m,i/2 ~ 700 GeV. The branching ratio for X2 - ► xh exhibits a broad peak of similar 
height to the lower values of tan (3, whereas the branching ratio for % 2 — > \Z is somewhat 
smaller. Benchmark I is located near the peaks of the X2 — > xh, X% branching ratios, and 
benchmark J is located in the region of mi/ 2 where the X2 & starts to rise. 

Panel (e) of Fig. 4 shows the branching ratios for tan (3 = 50 and fi > 0. Here we note the 
strong dominance of X2 — ► tt decays, the approximate constancies (at relatively low levels) 
of the branching ratios for % 2 — > xh, Z above their respective thresholds, and the strong 
suppression of % 2 — > ££ decays. Benchmark L has branching ratios that are typical for this 
value of tan (3. 

Finally, panel (f) of Fig. 4 shows the branching ratios for tan (3 = 10 when fi < 0. In this 
case, there is dip structure in the branching ratio for X2 which rises monotonically 

with TOi/ 2 , while the branching ratio for X2 ~^ rr is always large. In this case, the peak 
branching ratio for % 2 — > xh is below that for % 2 — > x%- 

The X2 branching ratios are rather different at the tip of the funnel for tan (3 = 50 - 
which has BR(%2 — ► xh) = 0.11 at the location of point M, by point K on the side of the 
rapid-annihilation funnel for tan/3 = 35 and \i < - where the mode X2 xh gives half 
of the total X2 decay rate, and at the focus points E and F - where the X2 decay rate is 
saturated by XQQ and by xh, respectively. 
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4 Production and Detectability along WMAP lines 



We now study the reaches at different accelerators in the number of observable super symmet- 
ric particles, and the changes in experimental signatures and topologies along the WMAP 
lines for different values of tan (3. 

4.1 LHC 

4.1.1 Sparticle Signatures 

In order to visualize the changes in the signatures of sparticle decays along the WMAP lines, 
we first compute the mean numbers of different particle species produced in sparticle decay 
chains at LHC. Large samples of inclusive supersymmetric events have been generated along 
the WMAP lines, using PYTHIA 6.215 [34] interfaced to ISASUGRA 7.67 [20] to compute 
the average numbers of different particle species produced per supersymmetric event, and 
the production cross sections of exclusive and inclusive sparticle production reactions. Fig. 5 
shows the average numbers of (a) Z, (b) h bosons, (c) r leptons and (d) trilepton events 
with three i (£ — e, /i) in inclusive sparticle decays at the LHC. The plots are all for /i > 
and display results for tan (3 = 10, 20, 35 and 50 along the corresponding WMAP lines. 

We see in panel (a) of Fig. 5 that the fractions of sparticle decays with a Z boson in the 
cascade is never large in the regions of the parameter space traversed by the WMAP lines, 
reaching about 0.1 Z°/sparticle decay at large and tan/5. 

The decays of sparticles into the lightest Higgs boson h, followed by h — > bb decay, provide 
another important signature at the LHC. These appear above the threshold for % 2 — > X^ 
decay at about m 1 / 2 =350 GeV, and may reach 7-10% per sparticle event at larger m^, as 
seen in panel (b) of Fig. 5. 

As seen in panel (c), the numbers of r leptons produced in sparticle events are always 
large, particularly at large tan/? where there are between one and two per event. Finally, we 
see in panel (d) of Fig. 5 that the fraction of trilepton events is always above 10%, and may 
attain ~ 40% at large m 1 /2 and small tan/3. 

These plots reinforce the importance of the r signature for sparticle detection at the 
LHC, within the CMSSM framework used here. Clearly the X2 — > xh, Z signatures are 
also interesting, but they may be quite challenging to exploit. However, we would like to 
emphasize that other decay patterns should not be neglected, and may be preferred in other 
supersymmetric models. For example, if the soft supersymmetry-breaking scalar masses for 
the Higgs multiplets are non- universal, X2 —> xh, Z decays may be more copious. 
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Figure 5: Supersymmetric event signatures at the LHC: the numbers of Z bosons (upper 
left), light Higgs bosons h (upper right), r leptons (lower left) and three-lepton final states 
per supersymmetric event are shown as functions of along the WMAP lines for four 
values of tan (3, all with fi > 0. These plots were obtained with PYTHIA 6.215 [34] interfaced 
to ISASUGRA 7.67 [20]. 
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4.1.2 Detectability at the LHC 

We now estimate the numbers of different species of supersymmetric particles that may be 
detectable at the LHC as functions of m^j along the WMAP lines for different values of 
tan /3, using the following criteria. 

• Higgs bosons: We generally follow the ATLAS and CMS studies of the numbers of 
observable bosons as functions of Ma and tan/3 [35]. In contrast to [3], here we also consider 
H/A-> X2X2 decays. 

• Gauginos: Our criteria for the observability of heavier neutralinos at the LHC have 
been refined compared to our previous publication [3]. In particular, we first compute the 
total numbers of \2 produced in all sparticle events and the branching fraction corresponding 
to a dilepton final state. We consider this to be observable at the LHC if the product is at 
least 0.01 pb, corresponding to 1000 events produced with 100 fb _1 of integrated luminosity. 
The lightest neutralino x is considered always to be detectable via the cascade decays of 
observed supersymmetric particles. 

• Gluinos: These are considered to be observable for masses below 2.5 TeV [36]. 

• Squarks: The spartners of the lighter quark flavours u, d, s, c are considered to be 
observable if rriq < 2.5 TeV [36], but we recall that it is not known how to distinguish the 
spartners of different light-quark flavours at the LHC. In general, we assume that the stops 
and sbottoms t, b are observable only if they weigh below 1 TeV, unless the gluino weighs 
< 2.5 TeV and the stop or sbottom can be produced in its two-body decays. 

• Charged Sleptons: These are considered to be observable when the mass splitting 
mi — m x > 30 GeV 11 and the inclusive production cross section, which includes their 
direct production and that in the decays of other supersymmetric particles, exceeds 0.1 pb, 
giving at least 10000 events with 100 fb _1 of integrated luminosity 12 . 

• Sneutrinos: These have not been considered as observable at the LHC, due to their 
large yield of invisible v v\ decays. 

We first note the key differences between the expected numbers of detectable MSSM 

particles at the updated benchmark points, shown in Fig. 3, and the analogous LHC analysis 

in [3]. Howver, we would like to caution the reader that it is impossible to be precie about the 

capabilities of the LHC (or other accelerator) without detailed simulations that go beyond 

the scope of this paper. 

11 This relatively conservative cut is to ensure that the decay lepton should be observable, and causes 
sleptons to be lost in some scenarios, particularly with the lower values of too at the updated benchmark 
points. 

12 This criterion is illustrative: in a more detailed study, one should also include backgrounds such as 
W + W~ , that may be very important. 
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• At point B, we now consider the heavier neutral Higgs bosons to be detectable via their 
decays into the xi and into r leptons, and the H ± — > ru decays should also be detectable. 

• We no longer consider the H ± to be observable at point C. 

• At point F, we now observe that the g has branching ratios of 16 % and 17 %, respec- 
tively, for decays into X2,3^, followed by branching ratios of 97 % and 99 % for \2 — > xh and 
X3 — > respectively. The g also has a branching ratio of 33 % for the decay into x ±: tb. 
We now consider that the X2,3 should be detectable at this point, but not the x ± - 

• At point H, supersymmetric particles are now well within the reach of the LHC, essen- 
tially because of the reduction in m^. We now consider the q, g, x and xi to be detectable 
at this point. 

• At point J, we now consider the X2 to be observable because of its large production 
rate, even though its decays generally include r leptons. 

• At point K, we now find that no sparticles are observable with our present criteria, 
although the lightest Higgs boson h is observable. 

A corollary of these changes is that supersymmetric particles appear to be detectable 
at the LHC at all the updated benchmark points except K and M, which are located in 
rapid-annihilation funnels. 

We now display in Fig. 6, as functions of mi/ 2 along the WMAP lines for different values 
of tan/3, the numbers of different types of MSSM particles that should be observable at the 
LHC with 100 fb _1 of integrated luminosity in each of ATLAS and CMS. The nominal lower 
bounds on mi/2 imposed by (dashed lines) and b — > (dot-dashed lines) are also shown. 
These each have some uncertainties, for example FeynHiggs has a quoted error of ~ 2 GeV 
in the calculation of rrih so that (specifically) benchmark B in the top right panel should not 
be regarded as excluded. 

We use the criteria explained above and generalize the results for the indicated benchmark 
scenarios that were shown in Fig. 3. In the top left panel, for tan/3 = 5, we note first that 
only one MSSM Higgs boson is expected to be visible. In the top right panel for tan (3 = 10, 
we see that the other Higgs bosons are also expected to be observable at low mi/2 in the 
neighbourhood of point B, for example via decays into X2X2- All of the Higgs bosons are 
expected to observable over the entire range of for tan/3 = 35,50 and \i > 0, but not 
for tan/3 = 10 and // < (bottom right panel). All of the charginos and neutralinos are 
expected to be observable at low m^, but only the lightest neutralino x at high m^. Since 
the sneutrinos decay invisibly and the rates for £ 2 production are inadequate, only the l\ 
are considered to be observable at the LHC for this value of tan/3, and only at low mi/ 2 . 
However, all the squarks and gluinos are expected to be observable anywhere along any of the 
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WMAP lines, except that for tan/3 = 50 and \x > 0, where sparticles would be unobservable 
at the end of the funnel, as exemplified by benchmark point M in Fig 3. 

In the cases of larger values of tan/3, shown in the other panels of Fig. 6, all the MSSM 
Higgs bosons are expected to be observable for all the WMAP ranges of m^. However, 
the observabilities of the different charginos and neutralinos resemble those for tan/3 = 10. 
Some of the £ 2 may be observable, at least for small mi/ 2 . As for tan/3 = 10, all the squarks 
and gluinos are expected to observable at the LHC anywhere along the parts of the WMAP 
lines that are shown, with the exception of the tip of the line for tan (3 = 50. 

We comment finally on the case of the metastable f\ that appears towards the end of the 
WMAP lines where AM = \m fl — m x \ < m T . There is no specific study of this scenario at the 
LHC, but it has been estimated in an analogous gauge-mediated supersymmetry-breaking 
scenario that the efficiency for observing a metastable particle with m = 640 GeV would 
be 25 % (increasing for larger masses), and that a signal should be detectable if the overall 
production cross-section exceeds 1 fb [37]. In view of the large sparticle production cross 
sections at the LHC, and the fact (see Fig. 5) that each sparticle event produces on average 
0.1 or more f\ particles per event towards the ends of the WMAP lines, we believe that a 
metastable f\ could be detected out to the ends of the coannihilation tails. On the other 
hand, it may not be easy to determine the mass accurately, as the peak measured from the 
relativistic 1/(3 factor broadens with increasing mass. For example, when m = 640 GeV the 
FWHM of the mass distribution is estimated to be about 250 GeV [37]. The error in the 
mass estimate may be reduced in a large statistical sample, but such an analysis would need 
to assume a realistic experimental environment in order to evaluate systematic errors. 

4.2 Detectability at e + e Linear Colliders 

Our criteria for the observability of supersymmetric particles at linear colliders are based on 
their pair-production cross sections. 

• Particles with cross sections in excess of 0.1 fb are considered as observable, thanks to 
their production in more than 100 events with an integrated luminosity of 1 ab -1 . 

• The lightest neutralino x is considered to be observable only through its production in 
the decays of heavier supersymmetric particles. 

• Sneutrinos are considered to be detectable when the sum of the branching fractions for 
decays which lead to clean experimental signatures, such as vt — ■> x ± ^ T — e ! r ) an d 
v T — > W + fi~, exceeds 15%. 

• The 77 collider option at a linear collider would allow one to produce heavy neutral 
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Higgs bosons via the s-channel processes 77 — * A and 77 — * H, extending the reach up to 
750 GeV for 0.5 TeV e ± beams and up to 1.5-2.0 TeV for 1.5 TeV e ± beams. A 77 collider 
may also be used to look for gluinos, but we do not include this possibility in our analysis. 

• Finally, we assume that a metastable fx could be detected at any linear e + e~ collider 
with more than 100 events, and note that the mass could be measured more accurately than 
at the LHC, by measuring the production threshold as well as 1/(3. 

We consider e + e~ collision energies y/s = 0.5 TeV, 1.0 TeV, 3 TeV and 5 TeV, and also 
the combined capabilities of the LHC and a 1-TeV linear collider. Comparing our present 
estimates of the physics reaches of linear e + e~ colliders of different energies with those in [3], 
we observe changes due both to the criteria adopted and to the mass spectra. We note briefly 
the principal changes. 

• The updated point F has no supersymmetric particles observable at 1 TeV centre-of- 
mass energy. This is in part due to differences in the ISA JET spectrum optimisation, where 
it now reproduces more closely that from SSARD. 

• On the other hand, the reductions in the parameters m^, m for point H make slepton- 
pair production possible already below 1 TeV. The same point now has one squark accessible 
at CLIC with centre-of-mass energy 3 TeV, and all of the squarks at 5 TeV. 

• Point M now has no squarks accessible even to CLIC operating with a centre-of-mass 
energy of 5 TeV. 

4.2.1 TeV-Class e+e Linear Colliders 

The capabilities of a linear e + e~ collider with y/s = 0.5 TeV are illustrated in Fig. 7 [21]. 
It can see some gauginos and sleptons if mi/ 2 ^ 500 GeV. Its capabilities are therefore well 
suited to the — 2- friendlier scenarios displayed in the left columns of Fig. 3. Moreover, 
the sparticles that it would see would complement those detectable at the LHC. Also, we 
recall that such a linear collider would be able to measure sparticle properties much more 
accurately than the LHC [3, 8]. 

A linear e + e~ collider able to deliver collisions at a centre of mass energy y/s = 1 TeV has 
the potential to complement significantly the LHC in studying the supersymmetric spectrum, 
as seen in Fig. 8. In particular, such a linear collider would observe most sleptons and 
gauginos as long as mi/ 2 is below about 600 GeV, as exemplified by benchmarks B, C, G, 
I and L. Also, it would detect at least one supersymmetric particle over the whole range 
of the WMAP lines, even up to their upper mi/ 2 limits exemplified by benchmark point H, 
except in the funnel cases: /i > 0,tan/5 = 50 and /i < 0,tan/3 = 35. The accuracy it would 
provide in the measurements of several sparticle masses would enable GUT mass relations to 
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be tested, thereby completing the exploitation of the LHC data. We note also that a 1-TeV 
linear e + e~ collider would be able to observe the lightest squark, the t±, at low values of 
mi/2- This possibility is exemplified by benchmark point B, as also seen in Fig. 3. 

The complementarity between the LHC and 1-TeV linear e + e~ collider is displayed clearly 
in Fig. 9, where we see that, together, they cover the majority of the MSSM spectrum over 
most of the TOi/ 2 ranges covered by the different WMAP lines. 

4.2.2 CLIC 

A 3-TeV lepton collider, such as CLIC, is expected to access almost all the sparticle spectrum 
for mi/2 < 700 GeV, as seen in Fig. 10. This would enable it, for example, to distinguish and 
provide detailed measurements of the different flavours of squarks that will not be possible at 
the LHC. Moreover, at larger mi/ 2 CLIC would be able to observe (almost) all the spectrum 
of Higgs bosons, sleptons, charginos and neutralinos [23] . CLIC will also be able to observe 
gluinos via squark decays at the focus-point benchmark E 13 . It would therefore provide full 
complementarity with the LHC along the full extension of the lines. 

Finally, we display in Fig. 11 the capabilities of a 5-TeV lepton collider such as CLIC. We 
see that all the MSSM particles are detected along all the WMAP lines, with the exception 
of the gluino that would generally have been seen at the LHC, and, for /i > and tan (3 = 50, 
squarks towards the tip of the corresponding WMAP line, including the point M. 

Despite the larger machine-induced backgrounds and beam energy spread, CLIC is ex- 
pected to perform measurements of the properties of accessible suspersymmetric particles 
with good accuracy [23]. Slepton and heavy Higgs boson masses can be determined to (9(1%) 
accuracy, also when accounting for realistic experimental conditions and resolutions. A sim- 
ilar accuracy can typically be obtained also for sparticles reconstructed through cascade 
decays, such as the the X2 discussed below. The availability of polarized beams is not only 
beneficial to increase the signal cross sections (such as in the case of slepton-pair production) 
but also as an analyzing tool. A combination of measurements of the stop-pair production 
with different polarisation states can be used to determine the stop mixing angle [23]. 

4.3 A Neutralino Case Study 

Many of the most interesting differences in the capabilities of the different colliders discussed 
above arise in the spectrum of charginos and neutralinos. In particular, while the LHC and 
a TeV-scale linear e + e~ collider are largely complementary, they may not be able to discover 
13 As already noted, we do not discuss here the possibility of observing gluinos in 77 collisions. 
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all the charginos, neutralinos and sleptons when mi/ 2 ^ 400 GeV, whereas CLIC can in 
principle observe all of them along all of the WMAP lines. To confirm these statements, we 
now compare the physics reaches of different accelerators for neutralinos in the (mi/2,m ) 
plane, with the conclusions shown in Fig. 12(a). 

For the LHC, we show Fig. 12(a) the region of the (m^, m ) plane for tan (3 = 10 in which 
dilepton structures due to the decay chains X2 -+££,£—> £x an d X2 — > XZ°, Z° —> £ + £~ are 
expected to be observable in the cascade decays of heavier sparticles. We consider £ — e, // 
only, though the case £ = r would also be interesting at large tan f3, and require ax dilepton 
branching ratio to exceed 0.01 pb, as discussed above. For our purposes, the relevant parts 
of the (toi/2,to ) plane are those allowed by the WMAP dark matter constraint. We see 
that here the LHC coverage extends up to mi/ 2 — 500 GeV, as also reflected in Fig. 6. 

In order to see how accurately the 3-TeV version of CLIC could measure heavier neutrali- 
nos, we have made a new simulation of the process e + e _ — > XX2 followed by the decay chains 
X2 -> U,i -> ix and X2 -> XZ°,Z° -> £ + £~, using SIMDET [38] to simulate the detector 
response and PYTHIA 6.215 [34] interfaced to ISASUGRA [20] to simulate the signal. Events 
with at least two leptons and significant missing energy were selected. Both slepton and 
Standard Model gauge-boson pair-production backgrounds were considered. Combinatorial 
backgrounds were subtracted by taking the difference of the pair e + e~ + events and 

the mixed e ± fi T events, and a sliding window has been used to search for a >5-a excess in 
the Mu mass distribution. The results shown in Figure 12(a) for tan/3 = 10 make manifest 
the extended reach provided by CLIC, covering all the range of mi/ 2 allowed by WMAP. 
At larger values of tan/3, decay chains involving the f become more significant, requiring a 
more detailed study that should include r reconstruction. 

To benchmark the 3-TeV CLIC capabilities for measuring the masses of heavy neutralinos, 
a representative point has been chosen at m-1/2 = 700 GeV, m = 150 GeV, /j, > and 
tan/3 = 10, along the corresponding WMAP line. Here, m X2 = 540 GeV, m x = 290 GeV 
and mi L = 490 GeV. As seen in Fig. 12(b), at CLIC the dilepton invariant mass distribution 
shows a clear upper edge at 120 GeV due to X2 — > £ + £~ l followed by £~l — > £~X, which can 
be very accurately measured with 1 ab _1 of data. However, in order to extract the mass of the 
X2 state, the masses of both the £l and x need to be known. In fact, making a two-parameter 
fit to the muon energy distribution, the masses of the JIl and x can be extracted with an 
accuracy of 3 % and 2.5 %, respectively, with 1 ab _1 of data [39]. Improved accuracy can be 
obtained with higher luminosity and using also a threshold energy scan. We therefore assume 
that these masses can be known to 1.7 % and 1.5 % respectively, which gives an uncertainty of 
8 GeV, or 1.6 %, on the X2 mass, accounting for correlations. This uncertainty is dominated 
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by that in the i and x masses. The accuracy on the determination of the end point alone 
would correspond to a xi mass determination to 0.3 GeV, fixing all other masses, which is 
not sensitive to the details of the beamstrahlung and accelerator-induced backgrounds. 

5 Conclusions 

We have seen in this paper how different colliders could provide complementary information 
about the MSSM spectrum in updated CMSSM benchmark scenarios and along lines in the 
(mi/ 2 ,m ) planes for different values of tan/5 and both signs of /i. In addition to direct and 
indirect laboratory constraints, we have implemented the constraints on cold dark matter 
imposed by WMAP and other astrophysical and cosmological data. 

We emphasize that the LSP is stable in any supersymmetric model in which R parity 
is conserved, such as the MSSM, in which case astrophysical and cosmological constraints 
on dark matter must be taken into account. In principle, the LSP might be some different 
sparticle, such as the gravitino or axino. Specific scenarios of this type include models 
with gauge-mediated supersymmetry breaking, and benchmarks for such models have been 
proposed elsewhere [8]. We recall, however, that the WMAP data on reionization of the 
early Universe disfavour models with warm dark matter [12], so that all such models must 
grapple with the constraint on cold dark matter provided by WMAP and other data. 

Alternatively, one might seek to avoid the cold dark matter constraints by postulating 
some amount of R violation. The collider signatures of any such model would differ from 
those of the MSSM if i?-violating couplings are not sufficiently small, and an additional set 
of astrophysical and cosmological constraints come into play if the i?-violating couplings are 
not very large. Any i?-violating scenario must address these issues. 

We have updated a set of thirteen CMSSM benchmark points proposed previously [3] in 
light of the current WMAP constraint on supersymmetric dark matter. Since this constraint 
has reduced the dimensionality of the CMSSM parameter space, we have also introduced 
'WMAP lines' in the (mi/ 2 ,m ) plane for different values of tan/3, and have studied the 
physics reach of different accelerators, by charting the decay signatures and numbers of 
observable particles. It is interesting to observe that some important signatures, such as the 
pattern of \2 decays, remain rather uniform along these lines. However, we note that different 
signatures may be recovered by relaxing some of the assumptions taken here. An interesting 
example is the appearence of large xi decay branching fractions into h° or Z° in models with 
non-universal Higgs masses. As these are of importance for defining the phenomenology 
of early runs at the LHC, we plan to consider them separately in conjunction with the 
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corresponding dark matter density constraints. 

The observation at the LHC of an excess of events containing large missing transverse 
energy could immediately indicate that supersymmetry is present, and that R parity is 
conserved, so that the Universe should contain relic LSPs. However, it is possible that the 
relic LSP density might not saturate the bound on cold dark matter from WMAP, since 
there could be other significant contributions to the cold dark matter, e.g., from axions or 
superheavy metastable particles. It is therefore interesting and important to study, within 
the CMSSM framework considered here, to what extent colliders can verify that the LSP 
'does its job' of providing the cold dark matter [13], yielding a relic density that is neither 
too small nor too large. 

The answer to this question requires a detailed study of the accuracy with which a given 
collider (or combination of colliders) can measure sparticle masses. General measurements 
of supersymmetric final states, based on the methods already established by ATLAS and 
CMS several years ago, involving, e.g., the end-points of kinematic distributions, should 
allow the masses of the visible sparticles to be measured with precisions of 10 % or better 
in many cases. This should in turn provide constraints on the fundamental parameters of 
upersymmetry at the level of a few percent in minimal models such as the CMSSM. At this 
point one should be able to test the compatibility of the region of parameter space favoured 
by the LHC with that preferred by cold dark matter experiments and hypotheses, e.g., with 
the WMAP lines. 

To go further, one needs to understand the sensitivity of calculations of Vl x h 2 to variations 
in the supersymmetric model parameters [32] . These sensitivities depend on their values, are 
different for different benchmark scenarios, and vary along the WMAP lines. Determining in 
general the accuracy with which fl x h 2 could be calculated would require detailed simulations 
going beyond the scope of this paper. However, existing studies of the previous version of 
benchmark point B yield a provisional answer in this case. 

The following sparticle masses are thought to be measurable with the indicated accuracies 
at the LHC [40] at SPS point la [8], which is equivalent for our purposes to benchmark B: 

m-g = 595.1(8.0), 

m~ qL = 540.3(8.8), m g - L = 520.4(11.8) , m~ bi = 491.9(7.5), m h = 524.5(7.9) 

m ih = 202.3(5.0), m ia = 143.1(4.8), m fl = 132.5(6.3), 

m x = 96.2(4.8), m X2 = 176.9(4.7), m X4 = 377.9(5.1) (13) 

Since these simulations were made using ISASUGRA, we fit the determinations (13) to the 
CMSSM parameters mi/ 2 ,m Q and tan/3. We assume that the sign of /i is known, and set 
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A = 0, recalling A Q does not in any case have a large impact on Q x h 2 , and find: 

m = (103 ± 8) GeV, my* = (240 ± 3) GeV, tan/3 = 10.8 ± 2. (14) 

We have adjusted the central value of m to correspond to the updated benchmark value 
for point B, and then propagated the above uncertainties into the calculation of Q x h 2 using 
Micromegas, with the result shown in Fig. 13. 

As seen in Fig. 13, the calculated relic density has a well-behaved distribution that is 
well fit by a Gaussian with 

Q x h 2 = 0.111 ±0.018 (15) 

We have made a similar analysis using SSARD, finding an identical central value of Q x h 2 and 
errors lo!o2o> -aooo an d -om7 associated with m , m±/* and tan/3, respectively. We conclude 
that collider measurements should, in principle enable the expected value of Q x h 2 to be 
calculated with interesting accuracy. Measurements at linear e + e~ colliders would increase 
the accuracy in (15) and also enable Q x h 2 to be calculated for other benchmarks where data 
from the LHC alone might be insufficient. 

The LHC and linear e + e~ collider measurements would also provide crucial input for 
predicting the cross sections expected in direct dark matter searches or the LSP annihilation 
rate in the Sun. For this, it would in general be important to measure the higgsino and 
gaugino components of the LSP, which may be possible at the LHC in some regions of the 
parameter space by measuring the decay modes of the heavier gauginos or by comparing the 
rates of x* — > ttx y i a a (virtual) slepton and x* — > -\h. 

Dark matter physics is another example of the complementarity of LHC and linear e + e~ 
colliders for exploring new physics in the TeV energy range. Very possibly, Nature does not 
choose the CMSSM, and it is desirable to formulate benchmarks for other possibilities [8]. 
However, the CMSSM remains one of the prime options for physics beyond the Standard 
Model that might be accessible to the next generation of colliders. We hope that this 
paper contributes to a better understanding of what they might achieve, given the present 
constraints on supersymmetry from laboratory experiments and cosmology, as updated here 
using information from WMAP. 
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Figure 6: Estimates of the numbers of MSSM particles that may be detectable at the LHC as 
functions of mu^ along the WMAP lines for /j, > and tan (3 = 5, 10, 20, 35 and 50, and for 
fi < and tan/3 = 10. The locations of updated benchmark points along these WMAP lines 
are indicated, as are the nominal lower bounds on mw 2 imposed by (dashed lines) and 
b — > S7 (dot-dashed lines). 
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Figure 7: Estimates of the numbers of MSSM particles that may be detectable at a 0.5-TeV 

linear e + e~ collider as functions of rax/2 along the WMAP lines for /i > and tan/3 = 

5,10,20,35 and 50, and for fi < and tan/3 = 10. The locations of updated benchmark 

points along these WMAP lines are indicated, as are the nominal lower bounds on 

imposed by m h (dashed lines) and b — > (dot-dashed lines). 
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Figure 8: Estimates of the numbers of MSSM particles that may be detectable at a 1-TeV 
linear e + e~ collider as functions of m\ji along the WMAP lines for > and tan/3 = 
5,10,20,35 and 50, and for fi < and tan/3 = 10. The locations of updated benchmark 
points along these WMAP lines are indicated, as are the nominal lower bounds on 
imposed by m h (dashed lines) and b — > (dot-dashed lines). 
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Figure 9: Estimates of the combined numbers of MSSM particles that may be detectable at 
a combination of the LHC and a 1-TeV linear e + e~ collider as functions o/m 1( / 2 along the 
WMAP lines for \i > and tan/? = 5, 10,20,35 and 50, and for \i < and tan/3 = 10. 
The locations of updated benchmark points along these WMAP lines are indicated, as are the 



nominal lower bounds on m 



i/2 imposed by (dashed lines) and b — > s'y (dot-dashed lines). 
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Figure 10: Estimates of the numbers of MSSM particles that may be detectable at the 3-TeV 

version of the CLIC linear e + e~ collider as functions of along the WMAP lines for 

fi > and tan/3 = 5,10,20,35 and 50, and for /i < and tan/? = 10. The locations of 

updated benchmark points along these WMAP lines are indicated, as are the nominal lower 

bounds on m 1 / 2 imposed by m h (dashed lines) and b — >• 57 (dot-dashed lines). 
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Figure 11: Estimates of the numbers of MSSM particles that may be detectable at the 5-TeV 

version of the CLIC linear e + e~ collider as functions of along the WMAP lines for 

fi > and tan/3 = 5,10,20,35 and 50, and for /i < and tan/? = 10. The locations of 

updated benchmark points along these WMAP lines are indicated, as are the nominal lower 

bounds on m 1 / 2 imposed by m h (dashed lines) and b — >• 57 (dot-dashed lines). 
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Figure 12: (a) Comparison of the sensitivities in the (mi/2,mo) plane of the LHC and CLIC 
at 3 TeVto X2 — ► X^ + ^ decay and (b) the dilepton mass spectrum from this decay observable 
at 3-TeV CLIC for the CMSSM point m 1/2 = 700 GeV, m = 150 GeV, y. > and tan (3 = 10 
discussed in the text. 
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Figure 13: The estimated accuracy with which Q x h 2 could be predicted on the basis of LHC 
data for the updated benchmark point B, using the projected experimental errors reported in 
[40], a fit to ISASUGRA 7.67 parameters and the Micromegas code. We find similar results 
using SSARD. 
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